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We present a field theory analysis of a model of two SU(2n)-invariant magnetic chains coupled 
by a generic interaction preserving time reversal and inversion symmetry. Contrary to the SU(2)- 
invariant case the zero-temperature phase diagram of such two-leg spin ladder does not contain 
topological phases. Only generalized Valence Bond Solid phases are stabilized when n > 1 with 
different wave vectors and ground-state degeneracies. In particular, we find a phase which is made 
of a cluster of 2n spins put in an SU(2n) singlet state. For n = 3, this cluster phase is relevant to 
^^®Yb ultracold atoms, with an emergent SU(6) symmetry, loaded in double-well optical lattice. 

PACS numbers: 75.10.Pq 


I. INTRODUCTION 

Interacting quantum systems with symmetries higher than the ubiquitous SU(2) one may display new interesting 
physicsii For instance, chiral spin liquids with non-Abelian statistics might emerge in quantum magnets with an 
extended SU(7V) symmetry!^ Although high symmetries are more difficult to maintain, there are several interesting 
possibilities of their realization. For instance, in their recent preprint Kugel et.al^ suggest condensed matter realiza¬ 
tions of the antiferromagnetic chains with high symmetry, such as SU(4) in the context of orbital degeneracyA.^ The 
second possibility is related to alkaline-earth or ytterbium ultracold atoms which have a peculiar energy spectrum. 
The ground state and a metastable first excited state have zero electronic angular momentum so that the nuclear 
spin I is almost decoupled from the electronic spin. This decoupling paves the way to the experimental realization of 
magnets with an emergent SU(A^) symmetry where = 2/ -|- 1 is the number of nuclear states which can be as large 
as 10 with strontium and ytterbium atomsi^Ai Recent experiments with such atoms loading into optical lattices have 
directly observed the existence of the SU(fV) symmetry and have determined the specific form of the interactions 
between atoms in the ground state and the excited stateJ^ii^ 

The simplest lattice SU(A^) model stems from alkaline-earth atoms in their ground state loaded into the lowest band 
of an ID optical lattice with a filling of one atom per site which best avoids three-body looses. For large repulsive 
interaction, the resulting SU(iV) symmetric magnet is described by the Hamiltonian: 

n = jJ2Pn,n+l, ( 1 ) 


where P„,„+i = <S')^5')^+i+const, with being generators of the su(A^) Lie algebra, is a permutation operator 

acting in the tensor product of the Af-dimensional Hilbert spaces n and n -\- 1. Model (HD is integrable by means 
of the Bethe-Ansatz approach.^^ For J > 0 its excitation spectrum is gapless and at small energies the dispersion 
is linear. Hence the theory is critical with A^ — 1 gapless relativistic modes; its universality class is the SU(A^)i 
Wess-Zumino-Novikov-Witten (WZNW) model with central charge c = N — 

When one considers chains coupled by generic interactions the integrability is lost though at the same time the 
physics becomes more interesting, as one can learn from the example of the problem of two coupled S = 1/2 
Heisenberg chainsi^^— The latter model has a phase diagram which includes topological phases containing seeds 
of superconductivityi^^— It also may used to describe experimentally observed confinement of fractional quantum 
number excitations existing for the single chain problem i^^i^^ 

In this paper we consider a one-dimensional version of the famous Kugel-Khomskii model^ namely a model of two 
weakly coupled SU(AI) chains (labeled 1,2): 


n = 




5 ( 1 . 1 ) 
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where are permutation operators acting on states on chain number a = 1, 2. In the simplest case the interaction 

V is just a permutation operator acting between sites of different chains. However, to select for our analysis any specific 
form of the interaction would be too restrictive. Instead we will use an alternative approach, namely, we will consider 
weakly coupled SU(A^) chains with general interaction, which may go beyond the two-spin exchange. As the result, 
the problem becomes then a perturbed conformal field theory (CFT) with Hamiltonian: 

n = W[SU{N)i]i + W[SUiN)i ]2 + ^ Aaf, y dx 0i,a(x)02.6(x), (3) 

a,b 

where W[SU{N)i]i ^2 corresponds to the WZNW Hamiltonian for SU(A^)i CFT on the corresponding chain and Oa,i 
and 02,& are operators acting on states of chains 1,2. In the weak coupling limit |Aab| << J one may take into 
account only relevant perturbations which greatly restricts the number of possible operators. On top of the SU(A^) 
symmetry, our choice will be further restricted by considering only (i) spatially homogeneous ladders and (ii) ladders 
with inversion symmetry. The latter consideration excludes operators with non-zero conformal spin. In this way 
a multidimensional phase diagram of the lattice model is projected on a manageable (in fact, in most cases just 
two-dimensional) phase diagram of the perturbed CFT. 

We will investigate the infrared (IR) physics of model (|31) when N = 2n. Our main conclusions are that the physics 
of the SU(2n) ladder with n > 1 in some crucial aspects is different from the SU(2) case. In particular, there are no 
topological non-degenerate phases. All phases contain local order parameters corresponding to Valence Bond Solids 
(VBS). There are crystals of two types: with wave vectors (tt/ti, 0) or (7r/n,7r) (2fc^ VBS) and (27r/n,0) {Akp VBS). 
The latter possibility did not exist for the spin S' = 1/2 ladder. It can be viewed as a cluster of 2n spins put in an 
SU(2n) singlet state. In the simplest case, i.e., N = i, this cluster corresponds to the plaquette phase found in the 
numerical analysis of Ref. of the SU(4) two-leg spin ladder with an antiferromagnetic interchain coupling. For 
N = 6, we find the emergence of a cluster phase of six spins, leading to trimerization which should occur in the phase 
diagram of the two-leg SU(6) spin ladder. The latter case is directly relevant to the insulating phase of double tube 
of ytterbium ^^^Yb ultracold atoms. 

The rest of the paper is organized as follows. In Sec. H, the continuum limit of weakly-coupled SU(2n) two-leg spin 
ladder is presented to identify the leading perturbation in Eq. The result is the continuum model ([H]). In Sec. 
HI this model is analyzed using the conformal embedding approach. As a result model ([5]) is expressed as a theory 
of Zjv parafermions coupled to SU(V )2 WZNW via some relevant operator (l27l) . This new formulation allows us to 
determine the nature of the phases. The resulting analysis is provided in Sec. IV and our concluding remarks are 
given in Sec. V. Finally, our paper is supplied with several appendices where some additional technical details are 
described. 


II. THE CONTINUUM LIMIT 

In this section, we determine the leading perturbation of model ([3]) by means of a continuous description of two-leg 
SU(2n) spin ladder with generic interactions. 

Let us start with the decoupling limit where the lattice model ([2]) reduces to two decoupled Sutherland models 
O- Its low-energy properties can be obtained by starting from the U(V) Hubbard model at 1/N filling with large 
repulsive U interaction . At energies below the charge gap ~ U, the SU(V) spin operators in the continuum 
limit are described byi ^^d^i^i 

Sf Jiiix) + Jiiix) + N^^\x) + e^^'^-^nfix) + .., (4) 

where 1 — 1,2 denotes the two decoupled chains and kp = ir/Nao (oq being the lattice spacing) since the underlying 
Hubbard model is 1/N filled (1 electron per site). In Eq. (|3]), ^ are the left and right currents which generate the 

SU(V)i CFT. They are defined in terms of the underlying left- and right moving Dirac fermions Lia,Ria as 

Jm = R\oT^pRw, jfL=LlT^pLip, (5) 

where a summation over repeated greek indices (SU(V) indices) is implied. In Eq. ([5]), is a generator of the 
su(V) Lie algebra in the fundamental representation, the 2kp and Akp parts of the spin density are related to 
primary fields of the SU(V)i WZNW model. They transform respectively in the V-dimensional fundamental and the 
N{N — 1)/2-dimensional antisymmetric representation of the SU(V) group. In particular, we have 

N/^ = X TrigiT^), 


( 6 ) 
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A being a constant, related to the charge degrees of freedom, which can be chosen real for a matter of convenience. In 
Eq. (in]), gi is the SU(fV)i primary field, or WZNW field, for the I th chain with scaling dimension {N — 1)/N. This 
operator transforms in the TV-dimensional fundamental representation of SU(A^) and can be expressed in terms of the 
fermionic operators through the non-Abelian bosonization approach d^i^^i^^^ 

gipc. - ( 7 ) 

being a bosonic field which captures the properties of the charge degrees of freedom of each chain 1 — 1,2. 

From the continuum representation ((4|), we derive a transformation of the SU(Af)i WZNW fields with respect to 
the one-step translation symmetry Tag : 

Tag ■ 51.2 —>■ ( 8 ) 

The next step of the approach is to find the leading perturbation in Eq. ©• One can identify it by means a of 
a symmetry analysis. To this end, let us recall what symmetry restrictions we adopt. First, we have the SU(A^) 
symmetry: gi —>■ UgiU\ with U belonging to SU(AT). Second, we consider the spatially uniform model. Then since 
matrix operator gi is not invariant with respect to translations ([ 8 ]) the perturbation can include only products of 515 ^ 
(or g'^ 52 )- From the inversion symmetry it follows that there are no operators with nonzero conformal spin such as, 
for example, Trg^dxg 2 - These considerations yield the following Hamiltonian: 


n = W[SU{Nh-,gi] + W[SU{N)i,g2] + Xi 


dx 


Tr(5i5^) +H.C. 


-I- A2 


dx 


TrgiTrgJ -H H.c. 


+ 


(9) 


where the dots stand for less relevant operators like marginal current-current interaction. For the case when the 
interchain interaction includes only two spins, i.e., Vn]r?^ = Jx '^a can check directly the form of Eq. 

(|9|). By substituting Eqs. (I4l6ll into the interaction term of Eq. ([Jj) and taking into account that the individual chains 
are described by the SU(A^)i WZNW model, we find model ® with A 2 = —Xi/N and Ai = J_lA^/2. It is important 
to note that the ratio of the naive continuum limit A 2 = —Xi/N is not universal and will be modified by higher-order 
contributions. For N = 2n matrix —g still belongs to the SU(A^) group. Then the transformation gi —gi leaves 
the WZNW part of (0) unchanged and changes the sign of the coupling constants of the perturbation. This fact will 
enable us to identify the corresponding parts of the phase diagram of model ([9]). 


III. CONFORMAL EMBEDDING APPROACH 

In the weak-coupling regime, two-leg SU(A^) spin ladders are thus described by a model of two coupled SU(A)i 
WZNW models perturbed by two strongly relevant perturbations with scaling dimension 2{N — 1)/N. The phase 
diagram results thus from the competition between these two terms. Though they have the same scaling dimension, 
the two perturbations are of very different nature. The one with coupling constant Ai is invariant under an SU(A)l 
X SU(A^)r symmetry: gi —>• Ui^giU^, I/l.r being two independent SU(A) matrices. In stark contrast, the second with 
coupling constant X 2 is only S\J{N) invariant. It turns out that at A 2 = 0, as it will be seen, theory Q is integrable 
and therefore it makes sense to consider this model at |Ai| >> IA 2 I and treat the A 2 -term as a perturbation. 


A. Fateev model: A 2 = 0 case 

In this respect, we now use the following conformal embedding which singles out the SU(A^) symmetry: 

SUiN)i X SU{N)i - SU{N)2 X Zn, (10) 

where the SU(A ^)2 CFT has central charge c = 2{N'^ — 1)/{N + 2) and Zjy is the parafermionic CFT with central charge 
c = 2{N — 1)/{N + 2 ) 4 ^ 1 ^ The latter captures the universal properties of the critical point of the Z^v generalization 
of two-dimensional (2D) Ising models. The low-temperature and high-temperature phases are described respectively 
by order and disorder parameters ak and Hk {k = 1,.., A — 1) which are dual to each other under the Kramers- 
Wannier (KW) transformation. The latter maps the Z^v symmetry, spontaneously broken in the low-temperature 
phase ((cT/c) 7 ^ 0 and (/ifc) = 0 ), onto a Z^v symmetry which is broken in the high-temperature phase where {g,k) 7 ^ 0 
and (cTfc) = 0. The order and disorder operators carry respectively a (k,0) and (0,fc) charges under the Zjy x Zjy 
symmetry: 


CTfe —7 under Zjy, Cfc —7 Cfc under Zjy 

g-k —7 under Z^v, gk ^ gk under Zjv, 


( 11 ) 








4 


with m = 0,...,-/V — 1. At the critical point, the theory is self-dual with a Zjv x symmetry and become 

primary fields with scaling dimension dk = k{N — k)/N{N + 2)i^ The Z^v CFT is generated by chiral right and 
left parafermionic currents 'i'kR,L l = ^ N-kR,L, k = 1,... ,N — 1) with scaling dimension = k{N — k)/N 

which are the generalization of the Majorana fermions of the Z 2 Ising model. Under the Zjv x Zjv symmetry, 
(respectively 'i’kR) carries a (fc, fc) (respectively {k,—k)) charge which means: 

'^kh.R —>■ under Zat 

'i>kL,R ^ under Zjv- (12) 


The next step of the approach is to observe that model (jH) at A 2 = 0, being SU(A)l x SU(A)r invariant, is 
independent of the SU(A )2 sector of the embedding (ITUl) but depends only on the Zjv parafermionic currents. One 
way to see that is to relate the Ai term of Eq. ®, i.e. Vi, in terms of the underlying Dirac fermions Ria, Lia of the 
continuum limit: 


Ui = - 


Ai 

Ai 

3? 


dx 


Aiq,I/2q 




H.c. 


= — TTT dx 


jtjR 


H.c. 


(13) 


where we have introduced a chiral SU(2)jv jL,R current. As shown in Ref. 1^ . there is a free-held representation of 
an SU( 2 ) 7 v current in terms of a bosonic held and the hrst Zjv parafermion current: 


3t = 

ZTT 

3l = ( 14 ) 


$-c 


($ic — ^ 2 c)lV^ being the relative charge 


bosonic held. For the right sector, we have a similar expression: 


Jr = 


Jr — 


-iy8lf7jv<i._cH^t 

27r 



(15) 


where the KW duality symmetry 'h iR — has been used for future convenience. 

We thus hnd that the Ai term in model is directly related to an integrable perturbation of Z^r parafermions 
introduced by Fateev with euclidean actio m^'^i^^ 


S = — A 


d^a;(4'iL^'ifi -f H.c .), 


(16) 


Szk being the action of the Z^r CFT and A = —AiA/47r^. This perturbation is invariant under the Zjv symmetry 
but explicitly breaks the Zjv symmetry as seen from Eq. m- 


B. SU(N)2 perturbed CFT 

Our next step is to express the A 2 -term in Eq. ® in the SU(fV )2 x Z^r basis. The expression of the two SU(A)i 
WZNW helds 31^2 in the SU(iV )2 x Z^r basis (N > 2) was obtained in Ref. [s^ We will justify their results from 
simple arguments based on symmetries. 

To perform such analysis we need the representation of the Z^v x Zjv symmetry in terms of the underlying Dirac 
fermions. This can be done thanks to the definitions ()14ll5p . Since 'I'il and 'I'ir have (1,1) and (1,-1) charges 
under the Z^r x Z;v symmetry, we find that the Zjv symmetry is implemented as follows on the fermions: 

Ela ^ ^ i?2a ^ e-*™/'^R2a, (17) 

while under Zjv we have: 

T2a ^ ^ R^c. ^ R 2 C.. (18) 










5 


From these results and the definition 0 , we deduce the transformation of the two original SU(-/V)i WZNW fields. 
Under the symmetry, we have 

whereas 51,2 are invariant under the Z^v symmetry. 

Now we are in a position to reproduce the results of Ref. [H which was based on CFT consistencies. First of 
all, since 51^2 transform in the fundamental representation of SU(iV), they should be related to the SU(iV )2 WZNW 
primary field G which transforms in the same representation and has scaling dimension {N'^ — l)/iV(A^ + 2) (see 
Appendix A). Since the scaling dimension of gi ^2 is 1 — 1/A^, we need an additional operator in the Zjy CFT with 
scaling dimension (iV— l)/A^(A^ + 2), i.e., cti, or the disorder fields /ri, /ij. One way to eliminate the ambiguity is to 
use the transformation of the different fields under the Z^v x Zjv symmetry. This suggests the following identification: 

{gi)a/3 Gapai 

{g2)aP ^ Gapul. ( 20 ) 

Indeed, the disorder operators /xi, g\ cannot appear in the decomposition since gi ^2 are singlets under the Zjv 
symmetry. The occurence of cti and a\ in Eq. dsni) are consistent with the transformation law of 51^2 under the Zjv 
symmetry m- We note that the results (1^ do not hold for TV = 2 which is a special case because the fundamental 
representation of SU(2) is self-conjugate. In that case, the expression can be obtained using the fact that SU(2)2 CFT 
is related to three decoupled 2D Ising models (see for instance Ref. [ 2 ^. One important consequence of the identity 
(|20| is that the one-step translation symmetry (|8]) becomes now: 

Tao : G ^ e^^’^/^G. (21) 

Finally, the identification (l20l) can be generalized for all SU(A^)i primary fields tpi which transform in the antisym¬ 
metric representation of SU(iV) described by a Young tableau with a single column and I boxes (1 < 1 < TV — l)i^ 
For the first SU(iV)i theory, these primary fields, i.e. (pu, are obtained by I fusion of gi by itself. Using the result 
(I 20 I) and the fusion rules of the Z;v parafermionic CFT, one can derive the correspondence between (pipi and SU(TV )2 
and Z^v primaries. If we denote <I>; the SU(A ^)2 primary field with scaling dimension 1{N + 1){N — l)/N{N + 2) (see 
Appendix A) which transforms in the 1th antisymmetric representation of SU(A^), we find: 

ipii ^ ^lai 

P 21 ^ ( 22 ) 

which is, of course, fine at the level of the scaling dimension since 1{N — l)/N = 1{N + 1)(TV — l)/N{N -1- 2) -1- 1{N — 
l)/N{N + 2). 

We are now ready to find to express the A 2 -perturbation of Eq. m in the new basis. Since cticti ~ cr 2 , we obtain: 

Tr giTr g+ Tr G Tr G+ : a 2 - Tr($adj)a 2 , (23) 

where $adj is the SU(iV )2 primary field with scaling dimension 2N/{N + 2), which transforms in the adjoint repre¬ 
sentation of SU(TV). In this derivation, we have used the definition of the adjoint primary fieldi^ 

Tr($adj) = Tt{G+T^GT'^) = i (^Tr G Tr G+ - ^Tr(GG+)^ , (24) 

and the identity 

— 2 — — SapSji^ . (25) 

As the result of Eq. (1^ we obtain the following expression for A 2 -perturbation, V 2 , of model (0: 

U 2 ^ A 2 J dx Tr($adj) (cr 2 + . (26) 

In summary, the low-energy effective field theory @ of weakly coupled SU(A^) Heisenberg chains can be reformulated 
in the basis with Hamiltonian: 

n = W[SU{N) 2 -,G] +nzY,-xJ dx + H.c.) + X 2 jdx Tr($adj) (^2 + 4) • 


(27) 
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The crucial difference between N = 2 and N > 2 cases is that for = 2 there is no (T 2 , crj operators and therefore two 
sectors of the theory decouple from each other. This does not happen for > 2 and this determines the difference 
in physics. Recall now that for N = 2n the spectrum of the original model ([9]) is invariant under a change of sign of 
the both coupling constants. Such sign change should be compensated by field transformations in Eq. (HZl), but we 
managed to find them only for N = 4p. 

IV. PHASES OF THE GENERALIZED TWO-LEG SU(2n) SPIN LADDER 

In this section, we will investigate the IR physics of model (|27)) to deduce the nature of the possible phases of 
generalized two-leg SU(2n) spin ladder. 


A. Field theory strong coupling approach 


To shed light on the possible phases, it might be interesting to first perform a strong-coupling approach directly to 
the continuum model (jHl) when both coupling constants Ai ^2 are of the order of the ultraviolet cut-off and |Ai| >> IA 2 I. 
In this respect, let us consider the euclidean action corresponding to Eq. 


S = S[SU{N)i-gi]+S[SU{N)i;g 2 ]+Xi 


dPx 


Tr( 5 ig^) +H.c. 


-I- A2 


cPx 


TrgiTrg^ -H H.c. 


(28) 


where the action of the S\J{N)k WZNW model is given by: 

S[SU{N)k;g] = d^x Tv (d^^g+d^g)+ T{g) 

r(ff) = ^ J^d^y {g+dagg+dfigg+d-fg), (29) 

g being an SU(A) matrix field and r(g) the WZNW topological term. 

The results of the strong-coupling approach depend on the sign of Ai and we assume N = 2n. 


1. Ai<0 

Then the minimization of the Ai term in action (l28l) gives gi = g 2 = G and the WZNW topological term ([29|) is 
doubled. The resulting effective action is therefore 

5eff = S[SU{N) 2 -, G] + 2 X 2 j d^x |Tr G\^. (30) 

Now, if A 2 < 0 then |Tr G\ should be maximal. For an SU(A) matrix it leads to the conclusion that G belongs to the 
center of SU(V), i.e. Zn, with G = J, fc = 0,1,..., V — 1 and I is the N x N identity matrix. The one-step 

translation Tag is spontaneously broken and the system has a finite order parameter with wave vector 2kF = T^/n. If 
X 2 > 0, as this would be the case for the two-spin interchain exchange ladder where A 2 = —Xi/N, we get the condition 
Tr G = 0. Here, as we will see, the order parameter has ikp wave vector. 


2. Ai > 0 


A similar approach leads to gi = —52 = G which is still an SU(V) matrix if N is even. The resulting model becomes 
then 


5eff = 5[.St/(V)2;G]-2A2 


d^x \TrGf. 


(31) 


For A 2 > 0 we have again G = which corresponds to 2kF order which is now staggered between the chains 

since gi = — 52 - For A 2 < 0 we get the condition Tr G = 0 corresponding to an Akp ordering. 
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B. Integrable point A 2 = 0. 


A weak-coupling approach can be performed by exploiting the fact that when A 2 = 0 model ([9]) becomes integrable 
and is related to the Fateev model dT51) . According to Ref. the integrable deformation of the Zjv parafermions 
m is a massive field theory with a mass gap A for any sign of its coupling constant A when N is even. The action 
explicitly breaks the Z^v symmetry and the exact spectrum consists of massive kink excitations that result from 
degenerate ground states labelled by an odd integer s = 1,3,... ,N + One can then average over high-energy 

degrees of freedom represented by the theory (IT 6 l) and obtain an effective field theory for the SU(A )2 sector in the 
low-energy limit E A. This theory describes magnetic excitations carrying quantum numbers of the SU(A) group. 
Using the result ([27|) of the previous section, we obtain the Hamiltonian density of this effective field theory: 

nes = 44 : + : lili :) + , (32) 

where 4 l chiral SU(A )2 currents and the mass M of the Z^r particle serves as the upper cut-off. The coupling 
constant of the adjoint primary field ^adj is 

r], = X2{a2+4)sM‘^^, (33) 


= 2N/{N + 2) being the scaling dimension of the adjoint primary field (see Appendix A). When A > 0, the vacuum 
average of the spin fields aj in the ground state s is known from Ref. EU: 


4j)s 

Qj = 


— (0s|(Tjj0s) 


sin 

7r(j + l)s 

N-\-2 

sin 

(wh) 


(M/4)^'‘^e^^ 


dt ( sinh(tj) sinh[(A — j)t] 
Jq t I sinh(At) sinh[(A-I- 2 )t] 


(34) 


where hj = j{N — j)/N{N + 2). In particular, from Eq. (IMl) we see that (cri)s = {aN-i)s = (cr^s- According to 
Eq. (1^ . this means that for Ai < 0, i.e., A ^ —Ai > 0, we have gi = g^, as we envisaged from the strong-coupling 
approach of the previous subsection. 

Let us now return to the problem of finite A 2 . Our task is to analyse the infrared physics of the SU(A )2 model (15^ 
perturbed by its adjoint primary field. Since it is a strongly relevant perturbation, we may expect that the magnetic 
SU(A) sector is always gapped. 

Fine details of the spectrum, however, depend on the sign and magnitude of coupling constant 77 s iSD. According 
to the result ((Ml) these depend on the ground state of model (ITBl) |0)s. We suggest that s is selected in such a way 
that the ground-state energy of (IMl) is the lowest. It is reasonable to think that this corresponds to largest spectral 
gaps. 

To determine a qualitative dependence of the spectrum on rjs it is useful to consider a direct semiclassical approach 
to the interacting Hamiltonian density (|M|) using the identity (IMl) : 


Hint = 77 «Tr G Tr G+, 


(35) 


where G is now an SU(A) matrix. When 77 s < 0, the minimization selects the center group of SIJ{N): 

G = exp(i27r/c/A)/, 77 s < 0, (36) 

where k = 0,..., N — 1 and the solution breaks spontaneously the one-step translation symmetry Tag (IMl) . The 
ground state is thus A-fold degenerate. 

When rjs > 0, the minimization gives then the condition that G is an SU(A) matrix with the constraint: Tr G = 0. 
The general solution for N = 2n takes then a Grassmanian form>i^ 

G = exp(i27rA:/A) [/■'"diag(l, ...1, —1,... — !)[/, 77 s > 0, (37) 

U being a general unitary U(A) matrix. As shown in Ref. |T^, within this semiclassical approach, model (IM)) with 
even N becomes equivalent to the Grassmanian sigma model on U{N)/\U{N/2) x U{N/2)] manifold with a trivial 
topological term 9 = 2^. This model describes a fully gapped phase. 

Although for both signs of rjs we obtain gapped spectra, the gaps for the Grassmanian sigma model are expected 
to be smaller. This becomes clear at large N » 1 since 1/A serves as a coupling constant for the sigma model and 
hence the gaps are exponentially small in A: m ^ exp(—constA). At the same time for rjs < 0 the gaps are algebraic 
in 1/A. 

From these results, obtained within the semiclassical analysis, we can now determine a phase diagram of the 
generalized two-leg SU(277) spin ladder. 








C. Ai < 0 


We first assume that Ai < 0 so that A ^ — Ai > 0 and we can use the result (l34)) for the vacuum expectation values 
of the CTj fields. As we have mentioned above, the coupling constant rjg (1551) of the low-energy theory (|55|) depends 
on the vacuum of model (1161) . The degeneracy is removed by a selection of r/g which yields the lowest energy ground 
state energy for model (1521) . 

To establish a relationship between s and rjg we recall that the vacuum expectation values (1341) enjoy the property: 
(o'j}s = (—= (o'])s since s is odd. Then substituting Eq. (1551) into Eq. (1551) we get for N = 2n: 


f}s=2k+l 


A 2 ^1 + 2 cos 


'7r(2k + 1) ■ 
- n -I-1 - 


). 


fc = 0 , 1 , ...n. 


(38) 


1 . A 2 < 0 

When A 2 < 0, the result (158)) suggests that the true ground state of model (1521) corresponds to s = 1. Indeed, 
(i) this state has the largest value of the coupling rjg, (ii) r]s=i < 0 which according to the semiclassical arguments 
presented above, corresponds to a larger gap than for a positive coupling. The semiclassical analysis conducted above 
predicts a spin gap phase with Wfold ground-state degeneracy which stems from the complete breaking the one-step 
translation symmetry Tag- Using the identification (1501) . we deduce the result: 

Trgi±Tr52-^TrG((Ti±aI), (39) 

and since ((Ti)s=i = (cr|)s=i ^ 0 from Eq. (1551) . we have: (Tr gi) = (Tr 52 ) ^ 0. A lattice order parameter can then 
be deduced by means of the result of Appendix B (see Eq. (1B71) with m = 1): 

{ 0:^1 in)) = e-*^ S^^SCn+i) - (Tr + Tr g,) ^ 0. (40) 

1=1 

Thus we have a uniform 2kp VBS phase which is in phase between the legs of the two-leg ladder and the wave vector 
is thus (tt/u, 0). This N = 2n-fold degenerate phase breaks spontaneously the Tag symmetry as expected from the 
semiclassical approach. 


2 . A 2 > 0 

In this case according to Eq. (1551) for even n the most negative coupling constant corresponds to a single vacuum 
k = n/2. For odd n there are two degenerate vacua k = {n± l)/2. 

Although the spectrum of model (1551) remains the same as for A 2 < 0, the order parameters change. Indeed, from 
Eq. (1551) we have 


and for fc = n/2 this vacuum expectation value vanishes. From the correspondence (1391) . we deduce that the 2kF 
order parameter vanishes as well: (Tr gi + Tr ( 72 ) = 0, in agreement with the strong-coupling approach. However, 
since = exp(i47rfc/A^)/ from Eq. (1561) and {cr 2 )s=n+i 7 ^ 0, for A^ = 4p we have instead the formation of a Akp VBS 
phase: 


2 

(OlkAn)) = SCnSt+i) - (Tr Tn + Tr :p 2 i) (Tr G^) ^ 0, (42) 

where we have exploited the result (lB7l) with m = 2 of Appendix B and Eq. (1221) . The Akp VBS phase has wave 
vector (27r/n,0) and a n-fold ground-state degeneracy. 

For n odd we have 




( 43 ) 
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and (cti) does not vanish for odd n. In this respect, one might expect the existence of a 2kp VBS phase described 
by the order parameter (liOl) . In contrast, the strong-coupling approach of Sec. IV A leads to the formation of a Akp 
VBS phase as in the n even case. One cannot exclude the occurence of a quantum phase transition between weak and 
strong coupling regimes when n is odd. However, the degeneracy between the s = n and s = n -I- 2 states might not 
be protected and with a result that the true ground state is a symmetric combination of the two states with opposite 
signs of (Ti so that the resulting average vanishes: (cri) = 0. Such scenario would support the Akp VBS phase as in 
the n even case. Numerical calculations on the two-leg SU(6) spin ladder are clearly called for to shed light on the 
issue of a quantum phase transition. 


D. Ai > 0 

We now turn to Ai > 0 case (A < 0). The "Ln model (fTBl) is still a massive field theory. Unfortunately, the vacuum 
expectation values of the order parameters (1341) for A < 0 are not known. However, when N = Ap (even n), one can 
perform the transformation '^il,r —>■ p to change the sign of A in model (ITSl) . From the fusion rules of the Zjv 

parafermionic theory^^ we have ui/ri ^ and ^ we deduce that 

CTl—:>iCTi, (72 —t — CT 2 , (44) 

and thus {(Ti)s = —(cr|)s. As shown in Appendix C, the Z 4 parafermions CFT can be described by a bosonic field 
theory with central charge c = 1. In particular, one can establish the mapping a 2 —> —<72 in the iV = 4 case directly 
by means of this bosonized description. 

The transformation ()44l) maps the spectrum of model (IS|) with couplings Ai > 0, A 2 onto the spectrum with couplings 
— Ai, —A 2 . The former one in Eq. (I44[) means that the order parameter in the region Ai > 0, A 2 > 0 is (Tr( 5 i — ^ 2 )) 
using the correspondence (IMl) . It leads to the formation of a staggered 2kp VBS phase with a iV-fold degeneracy: 

{Olffin)) = ^ (Trsi - Trg^) ^ 0. (45) 

;=i 

The emergence of a staggered phase is consistent with the strong-coupling approach in the Ai > 0, A 2 > 0 region. 


/ 

'^2 

0000 

f. . 

. .f. . . .). . 



♦ » ») 

f. . . .if. . . .i 

0000 




FIG. 1: Phase diagram of model m for A = 4 by means of a field-theory analysis when |Ai| >> IA 2 I and a mean-field approach 
for IA 2 I >> |Ai|. A similar phase diagram is obtained in the even n case where the plaquette phase is a n-site rectangnlar 
singlet phase. 

In contrast, the region Ai > 0, A 2 < 0 is occupied by Akp VBS with a V/2-fold degeneracy as in Eq. (11^ . Eor 
V = 4, it corresponds to a plaquette phase which is known to appear in the standard two-leg SU(4) spin ladder 
A summary of the results obtained for even n can be found in Fig. 1. 

When N = ApA-2, the vacuum expectation values (IMll of the Fateev model (HSl) with A < 0 are not known so we 
cannot conclude as before unfortunately. However, from Eq. ([5]), we saw that the change of sign Ai ^2 —>■ —Ai ,2 can 
be absorbed by the redefinition gi —>• —gi. The spectrum of model dH) for Ai > 0, A 2 is thus related to that of with 
couplings Ai < 0,—A 2 . We thus expect the formation of a staggered 2kp (respectively Akp) VBS phase when A 2 > 0 
(respectively A 2 < 0) as in the even n case. 
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E. Mean-field analysis in the IA2I >> |Ai| regime 

So far, we have exploited the existence of an integrable line A 2 = 0 in model dH) to find the possible phases of 
generalized two-leg SU(2n) spin ladders. As emphasized before, such approach is valid in the |Ai| >> IA 2 I regime. 
When Ai = 0, the strongly relevant perturbation in Eq. (I?7l) is not an integrable and is thus difficult to analyse. 
In the regime IA 2 I >> |Ai|, we expect the formation of spectral gaps in the problem. One can then investigate the 
resulting IR physics in this regime by means of a simple mean-field decoupling approach: T-L^nf = T~Li -I- 'H 2 with: 

"Hi = Uif, + X 2 J dx (Tr($adj)) ( 0-2 -f 0 - 2 ) 

^2 = N + 2^ ' : ) + A 2 J dx {a 2 + ctI) Tr(<i)adj). (46) 

The perturbation in the Zw sector with the a 2 spin field with scaling dimension 2{N — 2)/N{N -|- 2) is more relevant 
than the one in the Fateev model (fTC)) . In the regime IA 2 I >> |Ai|, one can thus ignore the latter perturbation. 
The field theory TLi is not integrable except in the N = 4 case where it is equivalent to the sine-Gordon model with 
/3^ = 27r/3 (see Appendix C). The 0-2 spin perturbation gives a mass gap to the Zn degrees of freedom and, when 
A 2 < 0, we have (ct 2 + CT 2 ) > 0 since (Tr($adj)) > 0. The perturbation in the SU(A^) sector, described by H 2 , becomes 
equivalent to model (I32p with rjs < 0. The one-step translation symmetry is then spontaneously broken and a 2]ip 
VBS phase emerges in the regime IA 2 I >> |Ai|. When Ai > 0 for instance, there is thus a competition between 
and 2kF VBS phases which is marked by the transition line in Fig. 1. The full determination of the transition lines 
is beyond the scope of our approach. 

A strong-coupling analysis of the lattice standard antiferromagnetic two-leg SU(2n) ladder is presented in Appendix 
D. In Fig. 1, it corresponds to the line A 2 = —Ai/iV with Ai —>■ 00 . For iV = 4 and N = 6, we find respectively plaquette 
and trimerized phases, i.e., 4:kp VBS phases. The lattice strong-coupling analysis confirms thus the prediction of our 
field-theory approach presented in Sec. IV when |Ai| >> IA 2 I. We thus conclude that the Akp VBS phase occurs 
in the phase diagram of the two-leg SU(4) and SU(6) spin ladders for a sufhdently strong J± > 0. Whether this 
phase is stabilized upon switching on a small interchain coupling or remplaced by a competing 2kp VBS phase is 
an interesting question which cannot be answered unfortunately within our approach. In this respect, large-scale 
numerical simulations are necessary to shed light on this intriguing possibility. When n > 4, an SU(2n)i quantum 
criticality is expected from the lattice strong-coupling approach of Appendix D. A quantum phase transition should 
occur between weak and strong coupling regimes of the standard SU(2n) spin ladder with J± > 0. 


V. CONCLUSIONS 

We have studied the zero-temperature phase diagram of the SU(2n) Sutherland antiferromagnetic ladder with weak 
but generic interchain interactions. We have found that for a translationally invariant ladder with inversion symmetry 
this phase diagram is controlled by only two parameters Ai ^2 in the continuum limit. These are coupling constants 
of the two most relevant operators allowed by the symmetries. Our field-theory analysis exploits the existence of an 
integrable massive field theory of Zjv parafermions along the A 2 = 0 line. In the vicinity of that line, i.e. |Ai| >> IA 2 I, 
our approach leads to the conclusion that the phase diagram contains only VBS phases characterized by local order 
parameters. In this respect, there are no topological phases as for n = I,i.e., a two-leg SU(2) spin ladder. There are 
two types of VBS phase with wave vectors (7r/n,0) or (7r/n,7r) {2kp VBS) and (27r/n, 0) {4:kp VBS). The Akp VBS 
can be viewed as a cluster of 2n spins put in an SU(2n) singlet with a n-fold ground-state degeneracy. For = 4, it 
corresponds to a plaquette phase which is known to appear in the standard two-leg SU(4) spin ladderi ^°d^ Our results 
for iV = 6 suggest that a cluster phase of six spins, leading to trimerization, should occur in the phase diagram of 
the two-leg SU(6) spin ladder. The latter case is directly relevant to the insulating phase of double tube of ytterbium 
^^^Yb ultracold atoms which can be engineered by considering double-well optical lattices4^ In this respect, it will be 
interesting to investigate numerically the latter system to confirm the existence of such a cluster phase for = 6 by 
means of numerical methods for SU(V) magnet as in Ref. 13. Finally, the field-theory analysis of two-leg spin ladder 
with SU(2n -|- I) spins is different and will be presented elsewhere. 
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Appendix A: Quadratic Casimir of SU(A) 

In this Appendix, we give some conformal data of some S\]{N)k primary fields which appear in the conformal embed¬ 
ding approach of Sec. III. The scaling dimension of a SU(iV)fe primary field which transforms in some representation 
R of the S\J{N) group is given byi^ 


Afl = 


2Cr 
N + k’ 


(Al) 


where Cr is the quadratic Casimir in the representation R. Its expression can be obtained from the general formula 
where R is written as a Young tableau: 


^ TL-tquj '^col 

Cr = = -[l{N- l/N) + E - E (A2) 

i=l i=l 

for Young tableau of I boxes consisting of rirow rows of length bi each and ricoi columns of length each. For 
instance, we get Cr = {N'^ — 1)/2N for the fundamental representation, Cr = N for the adjoint representation, 
= k{N + 1){N — k)/2N for the kth basic antisymmetric representation made of a Young tableau with a single 
column and k boxes, and Cr = N — 2/N + 1 for the symmetric representation with dimension N{N + l)/2. 


Appendix B: Lattice representation of SU(A)i primary fields 


We investigate here the lattice representation of SU(A)i primary fields which transform in the basic antisymmetric 
representations of SU(A^). The results are important for the identification of the phases of the generalized two-leg 
SU(2n) spin ladder. 

Let us consider an S\J(N) spin chain in the fundamental representation of SU(A^). As recalled in the introduction, 
the model is integrable and displays a quantum critical behavior in the SU(A)i universality class with central charge 
c = N — At low-energy, the lattice S\J{N) operators in the continuum limit are described by d^i^^ 


N-l 


— Jr + Jr + E 


^i2mkpx j\^A 


(Bi) 


where kp = n/Nao, x = nao, and Jr r are the left and right SU(A^)i currents. The 2mkF parts of this decomposition 
are related to the m = 1,..., A — 1 SU(A)i primary field with scaling dimension m{N — m)/N which transforms 
in the antisymmetric representation of SU(A) made of a Young tableau with a single column and m lines: 


Ni = Tr($™T;^), 


(B2) 


where are SU(A) generators in the mth basic antisymetric representation of the SU(A) group and Am non- 
universal real constants. It is interesting to get some lattice interpretation of Tr in terms of the original lattice 
SU(A) spin operators. In this respect, let us introduce lattice 2mkF bond SU(A) operators: 

= (B3) 

with m = l,...,A—1. The continuum description of these operators can be obtained from the identification (EJ: 

[Jl{x + ao) + Jr{x + Oq)] -f 6^^ [j^^) + Jr{x)] N^{x + Oq). (B4) 

At this point, it is useful to recall the defining operator product expansions (OPE) for the SU(A)i primary fields4S 

1 


(^) (w,w) - 2TTiz-w) 

Using this result together with the definition (IB2I) . we get the following OPE: 


(B5) 
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where Cm 
the lattice 


[Jl{z) + Jr{z)] N^{w, w) + N^{z, z) [J^(w) + Jr{w)\ 

■ ~ A r 1 1 ^ 


^mCm ( ^ 


27r 


z — w z — w 


Tr <^m{w,w), 


(B 6 ) 


is the quadratic Casimir of the mih antisymmetric representation of the SU(-/V) group. We then deduce 
representation of Tr in terms of the 2mkp bond SU(7V) operators: 

OlmkM) = ^ Tr <^m{x). (B7) 


Appendix C: A special case of the SU(4) group 

For TV = 4 one can exploit the fact that SU(4) group is isomorphic to 0(6). The Kac-Moody algebras SU(4)/c and 
0(6)/c are equivalent and for k = 1 and fc = 2 it is possible to employ Abelian bosonization and Majorana fermions 
techniques 


1. Z 4 Fateev model and sine-Gordon model 

In particular, the Z 4 parafermion CFT has central charge c = 1 and admits a bosonized description in terms of 
a compactified boson field y at radius R — y^3/27r (x ~ X + which lives on a Z 2 orbifold: x ~ ~X-“ Some 
primary fields of the Z 4 CFT have a simple bosonic representation in terms of a vertex operator as the 172 held and 
the first thermal operator ei: 


0-2 cos(\/27r/3 x), ei -- cos(y/ 87 r /3 x), (Cl) 

which have respectively scaling dimension 1/6 and 2/3 as it should. The spin field operators with scaling 

dimension 1/8 are related to the twist fields of the c = 1 Z 2 orbifold CFT4S 

Using this bosonization approach, the Fateev model (IT 6 I 1 can be shown to be equivalent to the = 67 r sine-Gordon 
theory with Hamiltonian density i^^i^^d^ 


■^24 = ^ {{dxxf + (dxX?) - Acos(v^x), 


(C2) 


X being the dual field of x- 

According to the exact solution of model (TIBI) , there are three degenerate vacua where CT 2 acquires expectation 
values (see Eq. (IMl) !: 


( 0 |cr 2 | 0 )^ - 1-f 2 cos( 7 rs/ 3 ), s = 1,3,5, (C3) 

corresponding to 2,—1,2. The sine-Gordon model (IC2I1 at /3^ = 67 r leads to the pinning of the bosonic field at the 
minima of the potential when A > 0 : (x) = P'\/27r/3, p being integer. Taking into account the radius of the boson 
which leads to the identification x X + we see that the sine-Gordon model (IC2() has three degenerate ground 
states with: (x) =0, i/27r/3, 2-\/27r/3. Using the bosonic representation (ICll) . we find that the vacuum expectation 
values of the a 2 operator give three values with two of them being equal as the exact result (|C3ll . Using this bosonic 
approach, we can also determine how result (IC3|) is modified when A < 0. The change of sign A —>■ — A can be obtained 
by the mapping: x X + 1^ + (?\/27r/3, q being integer. The value of q is fixed by the requirement that the first 

thermal operator should be invariant under the transformation since the Z 4 symmetry should not be broken when 
A < 0. The bosonization result (ED gives < 7=1 and therefore: 

X X + \/ 37 r/ 2 . (C4) 

We thus deduce that under the change of sign A —— A, the CT 2 operator transform as: tT 2 —>■ — a’ 2 - We recover the 
result (1441) from the bosonization approach for A = 4. 
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2. Abelian bosonization approach in the N — i case 


Using this bosonization description of Z 4 parafermions, one can investigate a free-held representation of model (l27l) 
when N = A. The central charge of the SU( 4)2 is c = 5 and this model might be described in terms of hve bosonic 
helds. 

First of all, the currents of the SU(4)i ^ SO( 6 )i CFT, which has central charge c = 3, can be expressed in terms 
of six Majorana fermions Xa- 

= ^XaXb, 1 < a < 6 < 6 . (C5) 

Then the SU( 4)2 currents being the sum of two SU(4)i ones, can be bosonized since two Majorana fermions can 
be expressed in terms of a single bosonic held: 


= iXaXb + irjarib = [cos(\/47ry:a) cos{-\/^(pb) + sin(\/4^(^a) sin(\/47ri^;,)] = 

TTOo 


HhUl cos[\/ 47 r((/?a - <Pb)], 

TTOo 


(C 6 ) 


where Kq are Klein factors. Since the center of mass of these bosonic helds drops out from this expression, we get 
the central charge 5. Thus we have a proof that SU( 4)2 can be bosonized. We suggest the following bosonic form of 
model (15^ : 



6 






(C7) 


The center of mass held 


= (C 8 ) 

^ a 

is decoupled from the interaction and is redundant. Here, however, it is important to obsverve that the helds cannot 
be the same as in Eq. (IC 6 I) because otherwise Tr$adj will not be local with respect to J. Using Eq. m and the 
identihcation (IC2p . we deduce a bosonic description of model (E7)) when = 4 with Lagrangian density: 


^ i - Ai cos(y67rx) - A 2 cos( V27r/3x) ^ cos 


a=l 


-i^a-^b) 


(C9) 


a>b 


From Eq. (IC3p . we see that the operator cos(-\/27r/3x) has three vacuum expectation values (two of which are equal 
to each other). The minimal ground-state energy is achieved when the sign of 7 = A2(cos(-\/27r/3x)) is positive. At 
I All >> IA 2 I we average over y held hrst and obtain the effective theory (jC7l) with 7 > 0. The latter can be studied 
either by 1/Wexpansion when we artihcially extend the summation over isotopic indices to N. Then we decouple the 
interaction via Hubbard-Stratonovich transformation and consider a saddle point: 


Ueff = ^ ^ + H.c) + 

' a 


(CIO) 


(7 ^ A 2 ). The saddle point represents iV-copies of the /3^ = 87 r /3 sine-Gordon theory; its spectrum consists of kink, 
antikink and one breather with mass of the kink’s mass. 


a. Confinement of the heavy partieles 

The structure of model (IC9I1 and especially the A 2 term critically affects the spectrum of the Z 4 Fateev model 
(which in the given case is equivalent to the sine-Gordon model with 0^ = 6tt (IC2I) '). Now we will try to discern the 
nature of the new excitations. 

As we have said, the operator cos(-\/27r/3x) has two degenerate vacuum values. When A 2 < 0 the values y = ±-\/ 7 r /6 
correspond to vacua with = 0. So one can suggest there are two types of kinks and antikinks of y-field: short 
kinks where V^^y shifts by 27r and long kinks where it shifts by 47 r. Both types of kinks are SU(4) neutral. 
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For A 2 > 0 it is not clear what the vacuum configuration is. It may be the vacua x = ±\/7r/6 and nonzero $0 or 
the single vacuum y = 0. Which situation is realized is determined by the energetics. In the first case there are short 
and long kinks as before, in the second one there is only one type of kinks, the one where \f^x shifts by Gtt. 

The second case corresponds to three particle confinement which has not been explored in the literature. This 
makes it worth commenting on which we do below. Averaging by small fluctuations of fields in (IC9I) we obtain 
the effective theory in the form of two-frequency sime-Gordon model: 

/Ipara = ^{^^l,xf “ ^1 cos(v^x) “ cos(v'27r/3x). (CII) 

At g = 0 the sine-Gordon model has only kinks and no breathers. When g > 0 the situation changes. The second 
term in Eq. (ICII|) has three vacua: Px = 0, 27r, —27r. In the presence of the last term the last two vacua become 
false. This leads to confinement of (0, 27r) and (0, —27r) kinks with formations of breathers (mesons). The kinks now 
are between 0 and Gtt vacua (hadrons) and are formed by confinement of three kinks (0, 27r), (27r, 47r), (47r, 6tt). Both 
mesons and hadrons are SU(4) singlets. Being topological excitations they will not experience decay from interaction 
with the SU(4) sector. 

The confinement of two kinks can be described in a standard way. The spectrum of mesons starts from above 2M 
threshould, where M is the kink’s mass in the /3^ = 6tt sine-Gordon model. 

The confinement of three kinks can be approximately described as follows. Only (0, 2tt) and (47r, 6tt) kinks interpo¬ 
late between the vacua with different value of (cos( y'^tt/Sx)). Kinks (27r,47r) do not change this expectation value. 
Therefore the linear potential exists only between the former kinks and the latter one can move freely in the space 
between the two. In the adiabatic approximation this movement contributes the amount 


6E = 


2M|a:i2P’ 

to the total energy. As a result the Schroedinger equation for the former kinks is 




4- = (A - 3M)4', 


(CI2) 


(C13) 


where M is the kink’s mass and rj ^ is the energy density of the false vacuum. 


b. The SU(4) sector 


Now let us return to the SU(4) sector. The kinks transform according to the fundamental representation of the 
SU(6) group, not the adjoint representation of the SU(4). Apparently, the above vector bosons are unstable and decay 
into the kinks, antikinks. 

How to get the SU(4) out of this? It is possible that since one has to project out the center of mass mode (IC8I) . 
the physical states are made of confined antisymmetric kink-antikink pairs which comprise 15-dimensional adjoint 
representation of the SU(4) group. 

The easiest task is to calculate the breather contribution to the {gg'^) since the neutral breather can be made as a 
linear conbination of breathers from different sine-Gordon copies. Such excitation just corresponds to the expansion 
of the cosines in (IG7I1 around their minima and hence does not contain the zero mode. So the contribution is 


Z 

_ p2 _ 3^2 ■ 


(CI4) 


Since the kinks and antikinks appear in pairs, the energy threshould for their emission is 2m or even higher if there 
is a real confinement and not just a condition of global neutrality. So it seems that the vector particles are indeed 
coherent as it has been found before. 


Appendix D: Strong-coupling analysis of the lattice two-leg SU(2n) spin ladder 

In this Appendix, we consider the standard two-leg SU(2n) ladder model ([2|) with V = (SU(2n) permutation 

operator between the two chains) in the strong-coupling regime when A —>■ oo. In that limit, the model is equivalent 
to a single SU(2n) spin chain where the spin operators belong to the antisymmetric representation defined by a Young 
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tableau with a single column of two boxes. For = 4, the physics of the latter model is well understood and a 
dimerized phase with a two-fold ground-state degeneracy occurs It corresponds to the plaquette phase of Fig. 1 
found in our field theory analysis in the weak-coupling regime. 

When 2n > 4, the nature of the phase might be inferred by means of a non-Abelian bosonization approach starting 
from the underlying U(2n) Hubbard model at a filling of 2 atoms per site {kp = 7r/nao)»^ At such a filling, there is 
an umklapp operator which couples the charge degrees of freedom with the non-Abelian ones: 

n 

Vu ll + H.C., (Dl) 

i=l 

where La^ and Ra^ are respectively the 2n left and right-moving Dirac fermions with ai = 1,..., 2n. Using the 
non-Abelian bosonization rule similar to Eq. we get a regular contribution: 

Vu ^ -gu cos ^V27rn$c) : (Trg)^ :, (D2) 

where is the charge bosonic field and g is the SU(2n)i WZNW primary field with scaling dimension 1 — l/2n. 
The operator : (Trg)^ : corresponds to the WZNW primary field with scaling dimension n/2 which transforms in 
the self-conjugate antisymmetric representation of the SU(2n) group with a Young tableau of a single column with n 
boxes. Assuming that we have a charge gap Ac (with a Luttinger parameter Kc < ^In), we have in the low-energy 
limit A <C Ac the effective interacting Hamiltonian for the non-Abelian degrees of freedom: 

-gu : (Trg)"^ :, (D3) 

which is a strongly relevant perturbation when n < 3, marginal for n = 4, and is irrelevant for n > 5. A straightforward 
semiclassical analysis gives rise, at ( 7 ci > 0 and n < 3, to an n-fold degenerate phase which breaks spontaneously the 
one-step translation symmetry: g = k = 0,..., n — 1. In this respect, a dimerized (respectively trimerized) 

phase is expected for n = 2 (respectively, n = 3), i.e., JV = 4 (respectively, N = 6). In particular, the trimerized 
phase for JV = 6 where the spins form a 6-site rectangular cluster, found in the weak-coupling regime, extends to 
the strong-coupling case. Interestingly enough, when n > 4, the interaction in Eq. (1031) becomes irrelevant and 
an SU(2n)i quantum criticality is then expected. For n > 4, we thus a expect a quantum phase transition, in the 
SU(2n)i universality class, for a finite value of the interchain coupling A. A recent variational Monte-carlo numerical 
analysis confirms these predictions and has shown that the SU(8) case displays a quantum critical behavior 
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